Interplane and intraplane heat transport in quasi two-dimensional nodal 

superconductors 
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We analyze the behavior of the thermal conductivity in quasi-two dimensional superconductors 
with line nodes. Motivated by measurements of the anisotropy between the interplane and intraplane 
thermal transport in Celrlns we show that a simple model of the open Fermi surface with vertical 
line nodes is insufficient to describe the data. We propose two possible extensions of the model 
taking into account a) additional modulation of the gap along the axial direction of the open Fermi 
surface; and b) dependence of the interplane tunneling on the direction of the in-plane momentum. 
We discuss the temperature dependence of the thermal conductivity anisotropy and its low T limit 
in these two models and compare the results with a model with a horizontal line of nodes ("hybrid 
gap"). We discuss possible relevance of each model for the symmetry of the order parameter in 
Celrlns, and suggest further experiments aimed at clarifying the shape of the superconducting gap. 

PACS numbers: 



I. INTRODUCTION 

Symmetry classification of possible gap struc- 
tures established the framework for separating con- 
ventional superconductors from their unconventional 
counterparts^^. Superconductors with the order pa- 
rameters that transform according to the trivial irre- 
ducible representation of the point symmetry group of 
the crystal are usually labeled conventional. If the order 
parameter transforms according to a non-trivial repre- 
sentation of the same group, the superconductor is called 
unconventional . 

In the latter class the symmetry of the order parame- 
ter is lower than the full group symmetry. This lowering 
of symmetry is usually believed to be related to strong 
repulsive Coulomb interactions and specific pairing mech- 
anisms (for example, due to magnetic fluctuations). As a 
result, determining the gap structure is one of the crucial 
steps in testing our understanding of the origin of super- 
conductivity in a given material. An important group 
of unconventional superconductors are those with order 
parameter that vanishes (has nodes) for some directions 
at the Fermi surface. 

Thermal conductivity is an exceptionally powerful 
probe for testing the shape of the gap in nodal supercon- 
ductors since only the unpaired electrons (with momenta 
close to the nodal directions) carry entropy. The presence 
or absence of the nodes, and sometimes their type (point 
or line) , can be inferred from the dependence of the ther- 
mal conductivity on the magnetic field and temperature. 
An additional test for the existence of line nodes is the 
so-called universality of the low-temperature coefficient 
of the thermal conductivity, ao = lim^^o k/T — The 
universal behavior refers to the relative insensitivity of 
ao to the concentration of the impurity atoms and to the 
details of the scattering on individual impurity centers^. 
Physically, since impurity scattering is pairbreaking, it 
generates near- nodal quasiparticles which can carry heat, 
and whose lifetime is, in turn, is limited by scattering on 



the same impurities. For linear nodes the two effects 
cancel nearly exactly. In the high-T c cuprates, for exam- 
ple, the quantitative agreement between theoretical esti- 
mates of ao and the measured thermal conductivity was 
established 8 -^. However, determining the location of the 
nodes in the momentum space from such measurements 
is much harder. 

The most direct tests measure the variation in the heat 
transport with the orientation of the applied magnetic 
field with respect to the nodal direction s 10 : 11 ' 12 : 13 ' 14 : 15 : 16 . 
Long before such experiments were attempted, it was 
proposed that the anisotropy of the thermal conduc- 
tivity along two different directions as a function of 
temperature allows to infer information about the gap 
structur o 17 ' 18 . Essentially, the measurement determines 
the predominant direction of the Fermi velocity for the 
nodal quasiparticles. Combined with the knowledge of 
the Fermi surface, measurements of the evolution of the 
anisotropy in the superconducting state impose stringent 
constraints on the possible loci of the nodes. 

This method was applied to UPty^, and the results 
argued convincingly for a line of nodes at k z = 0. The 
main reason this experiment did not uniquely determined 
the gap structure is that theoretically expected results 
depend on the detailed shape of the Fermi surface and 
the superconducting gap, i.e. on the basis functions for 
the particular representation^. This detail apart, it is 
believed that the measurement provides a good test for 
the horizontal vs. vertical linear nodes. 

Very recently the anisotropy of the thermal conductiv- 
ity along incquivalent crystalline directions was measured 
in heavy fermion Celrlns by Shakcripour et ali^I. The 
ratio of the c-axis to the in-plane thermal conductivity, 
n c /n a , rapidly decreases in the superconducting state. 
Similarity of the anisotropy evolution to that for UPt3, 
where equatorial line of nodes is believed to exist, led the 
authors of Ref. Hi] to suggest that the superconducting 
gap in Celrlns also has a horizontal line of nodes. This 
is surprising since CeCoIns, a close relative of Celrlns, 
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has vertical lines of node o 15 ' 22 , and the prevailing belief is 
that the two compounds arc quite similar (although some 
evidence points to different origin of superconductivity in 
the two systems 2 ^). 

Consequently, in this paper we re-visit the analysis of 
the temperature dependence and the anisotropy of the 
thermal conductivity in nodal superconductors. Since 
both de Haas - van Alphcn measurements, and band 
structure calculations show that Celrlns has an open, 
quasi-two dimensional piece of the Fermi surface with sig- 
nificant /-electron contribution 2 ^, our main focus here is 
on such anisotropic systems. We discuss the implications 
of the measurements for the symmetry of the order pa- 
rameter, and compare the anisotropy of the thermal con- 
ductivity for several models relevant not only to Celrln5, 
but also to other systems with a quasi-two dimensional 
Fermi surface. 

The remainder of the paper is organized as follows. In 
next section we introduce the basic experimental facts 
and theoretical considerations, formulating the models. 
The subsequent sections are devoted to the analysis of the 
proposed models. We then critically compare our results 
with experiment, and propose further measurements to 
test the gap symmetry in Celrlns. 



II. EXPERIMENTAL BACKGROUND 

Celrlns is a member of the so called 115 series, which 
also include CeCoIns and CeRhlns . While in the Rh sys- 
tem the /-electrons of Ce remain localized and undergo 
antifcrromagnctic ordering, both Ir and Co compounds 
are paramagnetic heavy fermion metals. In both systems 
de Haas-van Alphen measurements and the band struc- 
ture calculations indicate that a major sheet of the Fermi 
surface is quasi-two dimensional, although the energy dis- 
persion along the c-axis is substantial, as evidenced by a 
very moderate anisotropy (~ 2) between the out of plane 
and in-plane normal state transport coefficients. In con- 
trast to CeCoIns, which is a highly unconventional metal, 
likely close to a quantum critical point at finite magnetic 
fiel d 25 ' 26 , Celrlns near the superconducting phase is a 
good Fermi liquid, and therefore its excitations should 
be adequately described by the calculations in the Fermi- 
liquid framework. 

Both Celrlns and CeCoIns are ambient pressure super- 
conductors, and Knight shift measurements indicate sin- 
glet pairin g 27 ' 28 . Specific heat and in-plane thermal con- 
ductivity measurements^, as well as penetration depths 
and spin-lattice relaxation rat o 27 ' 28 , show the existence 
of line nodes. In CeCoIns there is overwhelming evi- 
dence for the vertical line nodes from the anisotropy of 
the thermal conductivity and specific heat under rotated 
magnetic fiel d 15 ' 22 , and from the tunneling into the An- 
dreev bound states^!. Given the strong similarities be- 
tween the Fermi surfaces it is tempting to conclude that 
the gap structure is similar in the two materials. There 
is, however, evidence pointing towards differences in the 



origin (and hence possibly in the type) of the supercon- 
ducting state of the Co and Ir compounds 2 ^. 

Authors of Refl2ll measured the temperature depen- 
dence of the thermal conductivity in Celrlns along two 
inequivalent crystalline directions, k c (T) (out of plane) 
and K a {T) (in plane), and made two significant observa- 
tions: a) The ratio R{T) = n c / n a is nearly temperature- 
independent above the superconducting transition tem- 
perature, T c , but is rapidly reduced with T at T < 0.5T C 
(Note that inelastic scattering yields a peak in thermal 
conductivity just below T c ~ 0.4 K, which may lead to 
the decrease of R(T) appearing very pronounced); b) The 
low T limit of the in-plane n a /T appears universal, but 
the interplane n c /T does not show the universal limit. 
Our goal below is to examine the possible origin and im- 
plications of this result. 



III. BASIC CONSIDERATIONS AND MODELS 

Independently of the choice of theoretical model, the 
results a) and b) above imply, prima facie, that the effect 
of superconductivity on the quasiparticles with the Fermi 
velocity predominantly in the plane is different from that 
for the quasiparticles moving along the c-axis: the states 
that carry entropy along the c-axis have a larger "effec- 
tive gap" than those carrying heat current in the plane. 
Several processes of different physical origin may lead to 
this behavior. In the following we do not consider inelas- 
tic scattering; this is appropriate at low temperatures, 
T ^_T C . The range of temperatures in Ref. [2l| where the 
inelastic scattering can be neglected is not immediately 
clear from the data, note, however, that, if the inelas- 
tic scattering is isotropic, our conclusions regarding the 
anisotropy of the thermal transport are not affected. 

Since the Fermi surface of Celrlns has several sheets, 
it is, of course, possible that the multiband effects are 
responsible for the observed behavior. However, without 
a microscopic theory detailing the gaps on each of the 
sheets, reaching reliable conclusions about the agreement 
between theory and experiment seems impossible, while 
phenomcnological multiband theory has too many fitting 
parameters to seriously constrain possible gap structures. 
We aim to construct a minimal model that captures the 
essential physics of the the experimental observations, 
and therefore restrict ourselves to considering one elec- 
tronic band. 

Since most band structure and dHvA analyses suggest 
an important role of the quasi-two dimensional (open 
along the c-axis) band, we choose such a band for our 
approach. We approximate it by a simple nearest neigh- 
bor tight binding expression 

k 2 + k 2 

e(k) = ^—^-2t(k)cos(k z c), (1) 
Ira 

where fcj are the components of the quasiparticlc momen- 
tum, fc, c is the lattice spacing in the z direction, m is the 
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effective mass, and t(k) is the intcrplanc tunneling ma- 
trix element. In all but one of our models we take t(fe) 
to be momentum-independent. We show below however, 
that the model with fc-dependent interplane tunneling 
that provides a viable path towards explaining the ex- 
perimental results, even with the assumption of vertical, 
rather than horizontal, nodes. 

We consider singlet order parameters, and begin with 
the simple model of vertical line nodes similar to the 
order parameter in CeCoIii5, introducing the azimuthal 
in-plane angle <p = arctaiiky/k x , 

• Model A: t(k) = t, and 

A(fe) = A o cos20. 

We show in the next section that this model is incom- 
patible with the experimental observations. Physically, 
the Fermi surface is cylindrically symmetric and the gap 
depends only on the azimuthal angle, <f>, while the z-axis 
component of the Fermi velocity is ^-independent. As we 
show below in Sec. [Vj these conditions ensure that the 
ratio n c /na = k zz /k xx , is temperature-independent even 
in the superconducting state. 

This constraint needs to be relaxed to explain the ex- 
perimental results, and we consider several possible mod- 
els. First, we show that for a tight-binding Fermi sur- 
face it can be expected, from microscopic considerations, 
that the gap also acquires a weak modulation along the 
c-dircction, while retaining vertical lines of nodes, and 
therefore consider 



• Model B 



t(k)-- 
A(fe) 



t, and 
= Ao(l- 



5 cos k z c) cos2(/>. 



While in this model there exists a temperature-dependent 
anisotropy, we find that this dependence is generally 
weak, and hence unlikely to provide a satisfactory ex- 
planation of the experimental observations (although we 
cannot exclude it). 

In the second approach to lifting the constraints of 
model A wc consider a situation with vertical line nodes, 
but where the nodal quasiparticles arc less efficient than 
the normal electrons (on average) in transporting heat 
along the c-axis compared to the xy plane. To model 
this we note that the measured and calculated Fermi sur- 
face of Celrlns is not rotationally symmetrioSi, and that 
the energy dispersion of the electrons along the z-axis 
clearly depends on the direction of the in-plane momen- 
tum. This dependence is expected in the D^h symmetry, 
and t(k) is determined by the structure of the overlap- 
ping wave functions for the quasi-two dimensional band. 
We therefore introduce the angle dependence into the in- 
terplane tunneling, preserving the crystal symmetry in 
the plane, and propose 

• Model C: t(k) = t + t\ cos 2 20, and 
A(fe) = A o cos20. 

For t\ > the interplane tunneling is smaller for the 
nodal quasiparticles than for antinodal ones, which en- 
sures the T-dependence of the anisotropy ratio. The sit- 
uation is somewhat reminiscent of that in the high-T c 



cupratcs where, in the absence of orthorhombic distor- 
tion, the interplane tunneling of the nodal quasiparticles 
is suppressed, i.e. t = 0—, and therefore the observed 
temperature dependence of the Josephson plasma reso- 
nance frequency^ 3 , and the c-axis penetration dept h 34 i 35 
is different from that suggested by the simple density 
of states power counting. We show that such a model 
gives a significant temperature-dependent anisotropy of 
the thermal conductivity. 

Finally, wc consider a model with horizontal, rather 
than vertical line nodes, similar to the hybrid E\ g gap 
proposed for UPt3, with broken time-reversal symmetry. 
The basis function of the representation is described as 
k z (k x +ik y ). For open Fermi surface we require periodic- 
ity along the ^-direction, and the appropriate basis func- 
tion is (k x +ik y ) sin k z c. The excitation spectrum and the 
thermal conductivity are only sensitive to the gap ampli- 
tude, |A(fe)| = k±\smk z c\, where k± = (fc 2 + fc 2 ) 1/2 . 
Since our consideration of Model B shows (see the anal- 
ysis below) that weak modulation does not change the 
anisotropy qualitatively, we ignore it, and consider 

• Model D: t{k) = t, and 

|A(fe)| = A |smfc z c|. 

We show that, this model also gives a substantially 
temperature-dependent ratio k zz /k X x- Moreover, wc 
point out that from the data of Ref . [2l| it is impossible to 
distinguish Model C from Model D, and suggest further 
measurements to probe the gap structure in Celrlns. 

We are now ready to consider each model in detail. 
In order to make connection with experiment, we largely 
focus on the thermal conductivity and its anisotropy 



R(T) 



n zz (T) 

l< XX 

(T)' 



(2) 



We begin by briefly looking at the normal state proper- 
ties. 



IV. NORMAL STATE 

The Fermi surfaces for the models A through D are 
shown in Fig. [T] In each case we transform to the cylin- 
drical coordinates and parameterize the Fermi surface by 
the azimuthal angle <j) = arctanfcy/fc x , and the c-axis 
quasimomcntum k z . The Jacobian of the transforma- 
tion from variables k x ,k y , k z to variables e, fc z , <fi is unity, 
and therefore no angle-dependence of the normal state 
density of states appears. We denote the Fermi energy 
Ep, and define the Fermi momentum at fc, = ±7r/2c as 
kp = \J2mEp. 

In models A, B, D above the tunneling t is <f>- 
independent. In that case, the radius of the Fermi surface 
is ko(k z ) — kpy/l + (2t/Ep) cos(fc z c), and the compo- 
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and therefore 



FIG. 1: Fermi surfaces for models considered in the paper. 
Left panel: Models A,B, and D. Right panel: Model C. Notice 
the similarity between the Fermi surface shown in that panel 
and the Fermi surface suggested for Celrlns in Ref. [24]. 



nents of the Fermi velocity are 
fco cos </> 



m 
fco cos 4> 



vpyl + {2t/Ep) cos(fc z c) cos 0,(3a) 



v F y/l + (2t/E F )cos(k z c) sin0,(3t>) 



v z = 2tcsm(k z c) — -jjj— (kpc)vp sin(fc z c), (3c) 

where vp = kp/m. In the normal state the anisotropy 
of the transport coefficients is given simply by 



Rn 



Ozz _ Kzz/T _ (V 2 Z ) 
&xx K xx/T ( v x ) 

where the average over the Fermi surface is 



(4) 



(A) 



d4> 

2^ 



" dxz 



A(<j),k z ), and Xz = k z c. (5) 



For the simple model above R n = {t/ Ep) 2 {kpc) 2 . 

We can obtain a very rough estimate of the parameter 
values relevant for Celrln 5 from the dHvA measurements 
and the band structure calculations^. For the quasi- 
cylindrical piece of the Fermi surface the typical radius 
is 1/5 of the Brillouin Zone size, kp ss 27r/(5a), where a 
is the in-plane lattice constant. In this material c/a ~ 
1.6, so that kpc ks 2. The normal state resistivity ratio, 
R n ~ 0.4, implying t/E F ~ 0.3. 

For model C the situation is slightly more complex. 
We still have 



v z = 2t{4>)csm{k z c) = 2i(0)csinxz, 



(6) 



tpti 



8E% 



(7) 



At arbitrary in-plane k the x-component of the velocity 
is 



fccos</> 2£i . , . , , 

v x = 1 — — sin <p sin 40 cos \z 

ui k 

and, consequently, at the Fermi surface, 



2 l l 
F 



sin 2 <h sin 2 ■ 



> cos 2 Xz 



l + 2t{<j})/E F cosXz 



(8) 



(9) 



Since we work in the regime t i/Ep <C 1, we can ignore 
the corrections from the term t{<p)/Ep in the denomina- 
tor as they appear only at the order (t/Ep) 4 , and write 



2 



1 



_3__tj_ 
16 E 2 F 



(10) 



For the conductivity anisotropy of 2.5 cited above, the 
maximal value oiti/Ep, when there is no to component, 
is approximately (ti/Ep) 2 w 4/15, which corresponds to 
the contribution of the second term of less than 5% to 
the in-plane transport in the normal state. Its effect in 
the superconducting state is even smaller, since the t\- 
dependent part of v x vanishes along the nodal directions. 
Consequently, we ignore the ii-dependent contribution 
to the in-plane transport hereafter. 



V. MODEL A: CONSTANT HOPPING, 
VERTICAL LINE NODES 

The thermal conductivity of a superconductor with the 
gap function whose average over the Fermi surface van- 
ishes is given b y 20 i 36 



~~T 



sech ( 



\2T I 



Ki(u,T), (11) 



Ki(u,T) 



1 



Re ( vf 



) 2 + |5| 2 -2|A k | 2 



U!lU>2 



Here Nq is the normal state density of states, and the 
renormalized frequency u~> = ui\ + iui-2 = & — where 
£ is the self-energy due to impurity scattering. For the 
model with Ak = Ao cos 2<p, and the hopping that is 
</>- independent, we immediately find for different orienta- 
tions of the heat current 



5E 



Ao 



2ru 

7T 

;i2 _ 



(12) 



21- 



K 



Ao 

UJ 



where E and K are the complete elliptic integrals, and 
Vx = vp, Tj z = (2tc) 2 . (13) 
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Therefore the anisotropy of the thermal conductivity in 
the superconducting state remains unchanged compared 
to the normal state, R s = R n = (t/ EpYikpc) 2 . More- 
over, in this model the universal (nearly independent on 
impurity concentration) low temperature limit, obtained 
from Eq.([12]), is- 



,4,5,6,20 



On , = hm — = , 

u ' 4 t^o T 6 



:E 



A 



VT 2 



(14) 



for both in-plane and out of plane directions. In Eq. (|14l) 7 
is the low-energy scattering rate, w(w = 0) = 17, and E is 
the complete elliptic integral. For clean superconductors, 
7 <C A , and we have 



1 



A 2 
^0 



-.E 



A„ 



A 2 



i +0 (k 



(15) 



so that a o i = r?j7riVo/(6Ao), universal and independent 
of the impurity scattering rate. 

This constant and temperature independent ratio 
R S (T) is in contradiction to the experimental measure- 
ments of Ref. and therefore model A does not provide 
a satisfactory description of the superconducting state of 
Celrln 5 . 



VI. MODEL B: MODULATED GAP 
A. Microscopic justification 

The symmetry classification does not uniquely deter- 
mine the functional form of the superconducting gap 
around the Fermi surface. The gap can be given by any 
combination of the basis functions transforming accord- 
ing to the chosen irreducible representation. For exam- 
ple, in a cubic system, both a constant gap and that 
varying as k x + k y + k\ have the same symmetry proper- 
tics. 

Of course, any fc-spacc variation of the gap beyond 
that imposed by symmetry requirements costs condensa- 
tion energy, and therefore, all other things being equal, 
the simplest basis functions often provide the most stable 
superconducting state. A conventional superconductor 
with a spherical Fermi surface strongly prefers a fully 
isotropic gap, while a strongly anisotropic quasi two- 
dimensional system with a circular Fermi surface and the 
dominant pairing in the B\ g channel is likely to have the 
gap with the k 2 — k 2 momentum dependence. This is the 
reason for using the familiar notation of s-wave for the 
former case, and d x ?_ y ? for the latter. 



On the other hand, superconductivity is an instabil- 
ity of the Fermi surface towards pairing, and therefore it 
is natural to expect that in real materials with complex 
Fermi surfaces the gap structure is also more complex. 
The correct gap can only be determined from the mi- 
croscopic theory of superconductivity. To solve Eliash- 
berg equations in superconductors with electron-phonon 
mediated pairing, Allen introduced the Fermi surface 
harmonics^!, which have the symmetry of the Fermi sur- 
face, are orthonormal, and are convenient basis functions 
for the expansion of the superconducting pairing field. 

In the absence of microscopic theory we use similar 
symmetry arguments, and go beyond lowest order basis 
functions. For simplicity, we continue to consider the 
d x 2_ y 2 gap symmetry, A(fe) oc cos 20; all our considera- 
tions apply equally well to other gaps with vertical line 
nodes, such as d xy . Choice of this symmetry determines 
the irreducible representation of the crystal point group. 
In a system with the Fermi surface open along the c-axis, 
we need to add to this point group the invariance with 
respect to translations by a reciprocal lattice vector along 
z. In the even (odd) channel the basis functions for trans- 
lation are cos nk z c (sinnfc z c), with n integer. Therefore 
quite 1 generally the singlet d. r ->_ ,,■> gap has the form 



A(fe) = A„ cos(nfc z c) cos 2</>. 



(16) 



Similar interaction for s-wave pairing was considered by 
Bulaevskii and Zyskir*^, and for d-wavc symmetry by 
Rajagopal and Jh a 39 ' 40 . Note that since even powers 
of k z , and therefore cos(nfc z c), transform according to a 
trivial representation of the tetragonal point group, the 
functions with different n belong to the same representa- 
tion and can be mixed. 

The degree of mixing is, of course, determined by the 
pairing interaction. We consider a separable model, 



V{k,k') = $(fc)$(fc'), 



where 



$(fc) 



V n cos(nk z c) cos 20. 



(17) 



(18) 



To illustrate the behavior of the model, and to be consis- 
tent with keeping only the nearest neighbor layer hopping 
in the energy dispersion, Eq.|T|), we truncate the expan- 
sion at n = 1, so that 



V(k, k') = Vq [1 + Ai(cosxz + cosx z ) + A 2 cos^z cosx 2 ] cos 20 cos 20'. 



(19) 
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To the same accuracy the gap is given by 
A(fc) = A [1 + 8 cos Xz] cos 20, 



(20) 



and is determined from the self-consistency equation (for 
pure superconductor in the weak-coupling limit) 



A(fc) = nT 



E 

|w„|<n 



dXz 
2tt 



V(k, fe')A(fe') 



2tt 



i + A3(fc') 

(21) 

Here V(k,k') = N V{k,k'), uj n = irT(2n + 1) are the 
fermionic Matsubara frequencies, and f2o is the cutoff 
energy for the pairing interaction. 

Formally, the gap given by Eq. (|2"D|) allows for a horizon- 
tal line of nodes when <5 > 1. This is, however, physically 
unlikely. In the interaction V(k, k') above, A2 is the rel- 
ative strength of the pairing in neighboring planes com- 
pared to the in-plane pairing, so that | A2 1 < 1- Generally 
we also expect Ai ~ t/Ep < 1— , so that the gap modula- 
tion along the c-axis is insufficient to produce horizontal 
lines of nodes. Inversely, if the interaction strengths in 
different channels in Eq.(|19p are comparable, we cannot 
truncate the expansion at the lowest order terms. There- 
fore, while in principle the symmetry considerations allow 
the gap in Eq. (|2"U|) to have both horizontal and vertical 
nodes, the model considered here can only be used if it 
yields small to moderate values of 5. 

Consider first the linearized equations for the tran- 
sition temperature, T c . Introducing dimcnsionlcss cou- 
pling constant g = NqVq, we find that T c and S satisfy 



1 + ¥ 



hi 



In 



7rT c 
2 lE Q a 



irT r 



(22a) 
(22b) 



Here je ~ 0.58 is the Euler's constant. Note that 
in the absence of the Ai term in the pairing interac- 
tion the equations for the temperatures of the transi- 
tion into the states Ao cos 20 and Ai cos Xz cos 20 decou- 
ple, and each is reduced to the standard BCS expression, 



T, 



M 



(27£j/7r)f2o exp(— 2/<7gjj), with the effective cou- 
pling constants in the cos nx z channels g and A2C//2 for 
n = and n = 1 respectively. For A2/2 < 1 this im- 



plies that Tc " 1 > Tc 1 ' , so that the simple cos 20 gives the 
dominant order parameter. 

For the generic case Ai ^ 0, however, the two channels 
are coupled, and the transition occurs directly into the 
c-axis modulated state A = Ao(l + Scosxz) cos 20. The 
transition temperature is 



2jE 



n exp(-2x /g), 



(23) 



where 



Xq 



(24) 



l + A 2 /2+DV2' 

D = (l-A 2 /2) 2 + 2A 2 . (25) 
We assumed that the coupling A2 is not too repulsive. 
The transition temperature depends on the magnitude, 
but not on the sign of Ai. On the other hand, the mod- 
ulation, 5 = ^(xq 1 — 1)/Ai, depends on the sign of Ai, 
and can be either positive or negative. 

The modulation amplitude, 5, is generally T- 
dependent. Consider the self-consistency equation at 
T<T C , 



l=gnT ^ / — / [l + A lC os Xz ] 



1 |<O 



\ui n \<n 



2 71 



2tt 



[1 + S cos x' z } cos 2 20' 
y/ivl + Al(l + 5cosx' z ) 2 cos 2 20' 

[1 + S cos x'z] cos 2 20' 
^JujI + A 2 (l + S cos x'z) 2 cos 2 20' 



(26a) 
(26b) 



In the linearized form considered above, all the terms in 
the integrand linear in cosx?, vanish by symmetry. On 
the other hand, below T CJ when A 2 term in the denomi- 
nator is important, these terms give a finite contribution 
to the gap equation. Moreover, in contrast to the con- 
ventional "one-channel" systems, the cutoff frequency fio 
cannot be removed from the self-consistency equations 
completely. As a consequence, S acquires temperature 
dependence. 



Of course, in the weak coupling limit this dependence 
is insignificant: in both integrals in Eg. (f2"S]) the greatest 
contribution is from the range Ao -C u>„ < flo, and there- 
fore essentially insensitive to the absence or presence of 
the gap. In most situations the dependence on the cutoff 
vanishes already for flo/T c > 10. This is seen from the 
sclf-consistcntly determined temperature evolution of the 
gaps in the two channels computed using Eqs. (|26[) and 
shown in the left panel of FigJ^l The corresponding plot 
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T/T T/T 

C C 

FIG. 2: Temperatures dependence of the gap parameters in 
Model B. Left panel: temperature evolution of the gaps in 
two channels for typical coupling values. Note the indepen- 
dence of the result on the cutoff frequency S7o: for different 
cutoffs Qo/T c = 10,100,1000 the curves lie on top of one 
another. Right panel: Temperature dependence of the gap 
anisotropy, 5 for different coupling constants. Note the tem- 
perature variation of the anisotropy in the regime of moder- 
ately large A2 = 1.6 and small Ai = 0.14, where we present 
the results for tt /T c = 10, 100, 1000. 



of the modulation, 8, in the right panel of Fig|5] demon- 
strates that the modulation is temperature-independent 
for most parameter values. 

We find that one exception is the somewhat artificial 
situation when a) the bare transition temperatures into 

the two states are very close, A2 ~ 2 so that « T"j , 
and b) Ai <C 1. While 8 is still temperature-independent 
for Q,q/T c 00, the temperature variations persist to 
large values of the cutoff, see FigfJl Part of the reason 
is that the value of 5(Ai — > 0, A2 — > 2) at T c obtained 
above depends on the order in which the limits are taken: 
setting A2 = 2 first gives 8 = \/2, while taking the limit 
Ai — > first for a finite value of 1 — A2/2 yields 5 = 0. 

Physically, at the point A2 = 2, Ai = the linearized 
equations, Eq. (|22[) . decouple, and give no information on 
the values of Ao and Ai = AqS. As is well known, the 
gap amplitudes are determined by the fourth order terms 
in the Ginzburg-Landau expansion of the free energy, or, 
equivalently, by the third order terms in the gap equa- 
tions. We find that inclusion of such terms gives three 
different possible solutions for Ao and Ai. Two of them 
are trivial, A 7^ 0, Ai = (lowest energy, as expected), 
and Ao = 0, Ai 7^ 0. However, we find that there is a 
third solution with A1/A0 = 5 = const, which has the 
highest energy of the three. At the same time for A2 = 2 
and any finite Ai 8 = const corresponds to the lowest free 
energy, which suggests a singular limit Ai — ► 0. This is 
confirmed by carrying our the expansion of 8 in 1 — T/T c , 
at finite Ai, where we find that the coefficients are sin- 
gular for A2 = 2. This suggests a strong dependence of 
the energies of the three stationary points on parameters, 
and on temperature. 
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FIG. 3: Gap amplitude as a function of the gap anisotropy 
for Model B at low temperature. The BCS value Aoo = 
2.14fcsr c . Note the smooth evolution of the gap amplitude 
across the value 8 = 1, where a horizontal line of nodes ap- 
pears. 



Away from the non-physical point A2 = 2, Ai = 0, 
the numerical solution of the self-consistency equations, 
Eqs. ptj)) shown in Fig fj] makes it clear that we can use 
the gap A = Arj[l + 8 cos \z] cos 2(f) with temperature 
independent S. However, at low temperatures Aq/T c is 
no longer at the BCS ratio of 2.14, but depends on the 
value of 8. This dependence is computed numerically, 
shown in FigjJH and is used in subsequent computation 
of the thermal conductivity. 



B. Thermal conductivity 

We use EcLlfTT)) with A = A [l + 8 cos Xz] cos 20 to 
evaluate the temperature dependence of the thermal 
conductivity in a quasi-two dimensional superconduc- 
tor. For our calculations we used the impurity scat- 
tering in the unitarity limit, with the normal state 
Tn = 0.007T C . While both k zz and k xx show the stan- 
dard dependence on the scattering rate for unconven- 
tional superconductors^, the anisotropy ratio is essen- 
tially insensitive to impurity concentration in the clean 
limit Tn -C T c . The results are shown in Figs. [4] and 
[5l One point of note is that the c-axis thermal conduc- 
tivity is even in 8, while the in-plane transport is sensi- 
tive to whether the largest or the smallest gap occurs for 
k z = 0, where the in-plane Fermi velocity is the greatest, 
see Eq.©. 

Of immediate interest to us is the temperature depen- 
dence of the anisotropy, R(T)/R n , which is shown in 
Fig. [6l The qualitative understanding of the tempera- 
ture dependence relies on the observation that the mod- 
ulation changes the gap most significantly in the regions 
Xz = 0, ±7t, where the quasiparticle velocity is strictly in 
the ab plane, with no component along the c-axis. Con- 
sequently, for 8 > the in-plane conductivity is reduced 
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compared to that of a superconductor with an unmodu- 
lated gap, while for <5 < it is increased. As a result, for 
5 > the anisotropy ratio k zz /k xx below T c is enhanced 
compared to the normal state value, while for 5 < it 
is reduced. Notice that for our model n zz is insensitive 
to the sign of d since points at the Fermi surface where 
cosx 2 = ±a (xz = oi, 7r — a) have identical values of the 
Fermi velocity, and therefore contribute equally to the 
c-axis thermal conductivity. 

The residual anisotropy in the T — > limit can be 
computed analytically once we recognize that, prior to 
integration over the z-component of the momentum, for 
each value of Xz the contribution of the kernel to the 
conductivity is universal in complete analogy to a two- 
dimensional d-wave superconductor with the gap A(xz)- 
Consequently, the ratio of the residual low temperature 
terms is (A = 2t/E F ) 
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1 + 5 cos Xz 
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;i + Xcos(x z ))dx z 
1 + 5 cos Xz 
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1 



(27) 



This result is shown in Fig. [5J The experimentally de- 
termined R s (T)/Rq is as low as 0.5 at temperatures of 
the order of 0.2T,£^. This can only be achieved for large 
negative values of the gap anisotropy, |<5| = — 5 > 0.8, 
when the gap nearly vanishes in the equatorial plane of 
the Fermi surface. This value is possible for sufficiently 
strong coupling Ai, but we consider it not very likely in 
Cclrlns. In this model both k zz and universal, 
albeit with different values, and therefore future doping 
studies of this material will be able to test the universal- 
ity of the z-axis thermal conductivity. 
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VII. MODEL C: 



0.0 



0.2 



0.4 0.6 

T/T 



0.8 



1.0 



FIG. 4: 

state. 



In-plane thermal conductivity for the modulated 



MODULATED INTERPLANE 
HOPPING 



A. Justification and general considerations 

In quasi-two dimensional systems the interplane ma- 
trix element depends on the overlap of the atomic or- 
bitals contributing to the bands close to the Fermi sur- 
face. Therefore, in general, the in-plane and out-of-planc 
motion of the quasiparticles are coupled, and, in writ- 
ing the tight-binding ansatz for the open Fermi surface 
along the c-axis, it is necessary to recognize that the ma- 
trix element t is a function of the in-plane direction of the 
momentum, <p. Moreover, simply looking at the quasi-2D 
sheet of the Fermi surface of Celrln.5 as obtained by band- 
structure calculations and by dHvA measurements^, it 
is clear that The Fermi surface is not rotationally in- 
variant, and therefore it is necessary to account for the 
angle-dcpendcncc of tunneling. For the Fermi surface of 



a general shape k 2 F = k 2 + k 2 — 2t cos Xz the difference 
between the maximal and minimal values of the radius, 
fcp(fc z ) is At. In the 115 series this difference is greater 
along the [110] direction than along the [100] direction, 
which suggest that t is 0-dependent 2 ^. 

In the normal state this dependence does little beyond 
replacing t with the appropriate average over the direc- 
tions, see Sec. IIV1 In a superconducting state with an 
anisotropic gap, however, the effect of the modulation of 
the intcrplane hopping can be much more pronounced. 
The best studied example is the high-T c cuprates, where, 
on the basis of band structure calculations, t(<f>) is ex- 
pected to vanish along the directions of the nodes in the 
gap22. Among the manifestations of this effect are the 
high power law in the temperature variation of the c-axis 
penetration depth (T 5 rather than T in the pure limit 
and weak dependence of the Joscphson plasma resonance 
frequency on temperature at T <C T c —. 

In the analysis of the thermal conductivity it is easy to 
see that any directional dependence t(4>) results in a devi- 
ation of the anisotropy ratio, R S (T), from its normal state 
value, R n , even for a superconducting gap that has pure 
d-wave form, A = Ao cos 20, with no fc 2 -dependcncc. 
This is simply because the quasiparticle velocity, v z , now 
acquires a dependence on the angle 0, which affects the 
evaluation of the thermal conductivity kernel, K z , see 
Eq. (|ll[) . Moreover, since the universal low temperature 
behavior is due to near-nodal quasiparticles, any sup- 
pression of the hopping matrix element in the vicinity 
of the nodes reduces the contribution of these quasipar- 
ticles to the c-axis transport, making it non-universal, 
while preserving the universality of the low temperature 
limit for the in-plane conductivity. Therefore assuming 
a modulated hopping may provide a route towards the 
explanation of the experimental results. 



where 



t(4>) =to + ti cos 2 20, 



(29) 



which is the simplest anisotropic form satisfying the 
tetragonal symmetry. This additional dependence does 
not affect the in-plane thermal conductivity, n XXl but 
does modify the out of plane k zz . As discussed in SeclIV) 
to a good approximation the normal state anisotropy ra- 
tio is 



-fin — — 



toti 3 ii 



El 



8E% 



(k F c) 



(30) 



For simplicity, and in agreement with the results of the 
previous section showing that small c-axis modulation 
of the gap does not yield significant corrections to the 
anisotropy ratio in the superconducting state, we con- 
sider a d-wave gap in this calculation. In agreement with 
the discussion above, we are interested in the situation 
when the intcrplane transport is suppressed along the 
nodal directions, and therefore take A = Ao cos 20. 

Neglecting the small corrections to the in-plane Fermi 
velocity due to the angle-dependence of t, as discussed in 
SccllVl we find 



B. Model and thermal conductivity 

Once again we adopt the model energy dispersion of 
the form 



e(k) 



2m 



(28) 
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(31) 

The kernel for the intcrplane conductivity is more com- 
plex. 
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where 



he{x) 
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he [x] 

hk(x) 
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3x 2 
K(x) - E{x) 



-E(x) 



3x 2 



-K(x), 



8a; 4 - 3a; 2 - 2 _, . 2(1 + a; 2 - 2x 4 ) „, . 

-E(x) + -A — J -K(x), 



15a; 4 
21 + x 2 



3 x 



—E(x) + 



2 + x' 
3a; 4 



15a; 4 
K(x). 



(34) 
(35) 
(36) 
(37) 



The T = limit, which is important for the comparison with the universal limit, can also be readily evaluated as a 
function of the low energy scattering rate, 7, by setting uj = 17. We reproduce the standard result for the in-plane 
conductivity, 



lim 

T^O T 



Kxx _ 7T N V 2 p 
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ey^TA 2 
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while for the intcrplane thermal conductivity we find 
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(38) 

(39) 
(40) 



In the clean limit, 7 <C Aq, the residual in-plane con- 
ductivity is universal, 



lim 1^. ~ ^ -^ oi; ^ 



(41) 



6 A ' 

while the c-axis conductivity, for tn/*i > 7/A0, is given 
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FIG. 7: Temperature dependence of the in-plane and the in- 
terplane thermal conductivity for the model with modulated 
hopping. The impurities are assumed to be in the unitarity 
limit, with V the normal state scattering rate. 
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which implies that the residual anisotropy ratio is 



Rq 
Rn 
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< 1. 



(42) 



(43) 



Even for a very moderate anisotropy ratio, t\/to = 2 
(which, for samples with 7/A0 ~ 0.1 is well within the 
range of the approximation), the anisotropy ratio at low 
temperature is already low at about 22%. 

The reduction of the anisotropy is even more pro- 
nounced for the cases of extreme anisotropy, to/ti <C 
7/A, when 



Rq 
Rn 



8 7" 
3A 2 



1 



to, 4A 
— In 

*i e-l 



< 1. 



(44) 



Consequently in this model the residual anisotropy ra- 
tio can be arbitrarily small depending on the purity. If 
there are no symmetry requirements for the hopping ma- 
trix element to vanish for the nodal directions, in most 
situations we do not expect to and t\ to differ by an order 
of magnitude; this suggests that the residual linear term 
is visible in the c-axis thermal conductivity, but may be 
significantly smaller than its in-plane counterpart. While 
the results above give a rough estimate of the effect of im- 
purities on the residual anisotropy, it is worth noting that 
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FIG. 8: Temperature dependence of the anisotropy of the 
interplane to the intraplane thermal conductivity. Note that 
the residual value is increasing with increasing to/ti. 



FIG. 9: Temperature dependence of the in-plane and the in- 
terplane thermal conductivity for the hybrid gap. The im- 
purities are assumed to be in the unitarity limit, with F the 
normal state scattering rate. 



disorder is likely to make the interplane hopping more 
isotropic, and therefore tends to restore the anisotropy 
to the the normal state value faster that what Eq. 
indicates. 



VIII. MODEL D: HYBRID GAP 



The temperature dependence of k xx and k zz is shown 
in Fig. IVIIIl and the corresponding anisotropy ratio for 
the hybrid gap is given in Fig. IVIIIl for the hybrid gap 
for different values of the impurity scattering. At low T 
the in-plane thermal conductivity is universal, while the 
interplane conductivity, 



So far we looked at the models with vertical line nodes, 
motivated by possible similarities between Celrlns and 
CeCoIn5, and attempted to reconcile them with the ex- 
perimental measurements. Of course, in the absence of 
any information about the nodal structure, a natural ex- 
planation for the observed anisotropy is that, in analogy 
to UPt3, the system has a horizontal line of nodes, and 
therefore the nodal quasiparticles do not contribute to 
the c-axis transport. We consider this model now. 

The hybrid gap belongs to the representation that 
transforms as k z (k x ± ik y ). The basis function for 
this representation over an open Fermi surface is (k x ± 
iky) sin^z, and therefore, for, if we take the weakly mod- 
ulated limit, kp x + kp y «const, the gap amplitude is 
|A(Xz, (j>)\ = Ao| sinxzl- If we take into account the four- 
fold, rather than cylindrical, shape of the Fermi surface, 
there may be additional small modulation of this gap 
with the component of the Fermi momentum in the xy 
plane as a function of z: based on our results for Model 
B in Sec lVIl we ignore these. Then the density of states, 
and the in-plane thermal conductivity for this model are 
identical to that for a system with vertical line nodes, 
while the interplane thermal conductivity kernel is given 
by 
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gives the residual anisotropy 
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This result is clearly non-universal. Notice that the de- 
pendence on the impurity scattering here differs from the 
7/A result obtained for UPt 3 in Ref . I20L The reason for 
the discrepancy is that for the hybrid gap and a Fermi 
surface closed along the c-axis the quasiparticles con- 
tributing the most to the c-axis conductivity are near the 
north and south pole. For the open Fermi surface these 
quasiparticles are absent, and the c-axis conductivity is 
reduced by an additional factor of 7/A0. It is also obvi- 
ous, from comparing this with the results of the previous 
section, that the residual anisotropy can be close for the 
models C and D, and may not provide the unequivocal 
distinction between the two. 
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FIG. 10: Temperature dependence of the anisotropy of the 
interplane to the intraplane thermal conductivity for the hy- 
brid gap. For comparison, we include results with deviations 
from the unitarity scattering limit. 

IX. DISCUSSION. 

It follows from experiment that the low-T value of the 
anisotropy of the thermal conductivity is below 25% of 
the normal value. The analysis above shows that model 
A (simple <i-wave with constant interplane coupling) does 
not give results that can explain the observed anisotropy. 
Model B (weakly modulated, along the c-axis, amplitude 
of the d-wave gap) shows some tendency towards the 
temperature-dependent anisotropy, but, in our view, re- 
quires fine tuning if it were to explain the measurements 
of Ref. [2TI- On the other hand, both model C (modu- 
lated interlayer hopping, vertical line nodes) and model 
D (horizontal line nodes) give very similar behavior of 
the thermal conductivity as a function of temperature, 
with a low residual value of the anisotropy k c / n a , and 
therefore may be relevant to Celrln 5 . Distinguishing be- 
tween the two based purely on the thermal conductivity 
data may not be easy, as is seen from the comparison in 
FigJTTJ which shows that, for sufficiently high anisotropy 
of ti/to, there is essentially no difference in the behavior 
of the ratio n c / n a as a function of temperature between 
the two models. 

There are notable differences between the results pre- 
sented for either model and the experimental observa- 
tions of Ref. |2j In experiment, the anisotropy ratio does 
not decrease below 0.5 for temperatures as low as 0.1T C , 
and there is no sign of saturation of the ratio below 0.1 5T C 
as in theoretical analysis. This is due to strong inelas- 
tic scattering that results in a peak in n c / n a at about 
0.3T C , but which is not included in our analysis here. 
Correspondingly, the low-T fit to experiment is done just 
below the peak, and the rapid decrease of the anisotropy 
may be related to the inelastic scattering. 

The low residual value of the thermal conductivity 
anisotropy implies, in model C, highly anisotropic tun- 



FIG. 11: Comparison of the temperature dependence for the 
hybrid and the d-wave gap. 

neling, t\ > 2to. Published band structure data do not 
allow to reliably extract this ratio from fitting the quasi- 
two dimensional sheet of the Fermi surface, but this is the 
task that perhaps should be attempted in near future. In 
the meantime, we propose several experiments that have 
the potential to provide additional information to resolve 
the question of vertical vs. horizontal line nodes. 

1. Direct measurement of the specific heat and/or 
thermal conductivity in the vortex state under a ro- 
tated applied field. This is, in our view, the best test 
for the existence of the line nodes. If, for the field 
rotated in the xy plane, no difference in the specific 
heat is found across the H-T phase diagram for the 
field direction along [110] vs. [100] direction, it is 
likely that there no vertical lines of nodes. If, on 
the other hand, such anisotropy is found and devel- 
ops in agreement with theoretical prediction a 41 ' 42 , 
this would be a strong evidence for c?-wavc, rather 
than hybrid, gap. 

At the same time, it is possible, although, in our 
view, unlikely, that there is gap modulation both 
in the plane and along the c-axis. Therefore mea- 
surements of the specific heat and thermal conduc- 
tivity when the field is rotated in the zy or xz 
plane should be used to eliminate this possibility, 
and directly probe for horizontal line nodes. Low 
superconducting transition temperature of Celrln5 
and inelastic scattering require that these measure- 
ments be done at temperatures below ~150mK. 

2. Measurement of the c-axis thermal conductivity in 
CeCoIn$. In the Co compound there is a gen- 
eral agreement that there are vertical line nodes. 
Measurements of the temperature evolution of the 
anisotropy n c / n a in that compound, and compar- 
ison with Celrln.5, while short of proof, would 
provide a test for the connection between the 
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anisotropy and the location of the nodes. Note that 
there is some evidence for remaining unpaired elec- 
trons in the superconducting state of CeCoIns (at 
least upon La-doping)^, which leads to a resid- 
ual metallic linear term for both directions of the 
heat current. However, as it was suggested that the 
fraction of these quasiparticles can be determined 
experimentally^, their contribution can be sub- 
tracted, and the remaining anisotropy used to test 
the differences and similarities between CeCoIn 5 
and Celrln 5 . 

3. c-axis penetration depth measurements. In anal- 
ogy to cuprates, if the interplane transport is sup- 
pressed along the nodal directions, the c-axis pen- 
etration depth should cither be not linear in T at 
T«T C (for t\ 3> to) or have a small range of linear 
T behavior, with the coefficient much smaller than 
that inferred from the density of states varying as 
E/Aq. In CeCoIns the c-axis penetration depth is 
linear in T, and observation of the non-linear be- 
havior in Celrln5 would point towards differences 
between the two systems. Unfortunately, we expect 
that, for reasons similar to those outlines for the 
thermal conductivity, it would be difficult to dis- 
tinguish model C from model D on the basis of this 
measurement. Negative result, however (linear, in 
T, A C (T)) would be difficult to reconcile with the 
thermal conductivity measurements. 

4. Doping studies of residual c-axis thermal conductiv- 
ity in Celrln§ . In Fig |12l we show the evolution of 
the residual thermal conductivity along the c-axis, 
in units of the expected universal value for that di- 
rection, as a function of the low-energy scattering 
rate, 7. Notice that the increase in k zz /T for the 
hybrid gap is slow, while for the modulated inter- 
plane hopping of model C it is faster. In reality 
we expect that addition of impurities will tend to 
make the interplane hopping more isotropic in addi- 
tion to introducing pairbrcaking, and therefore the 
increase in the residual value of k zz /T will be even 
faster than that suggested by Fig[T2J We therefore 
predict that, if model C of vertical line nodes with 
modulated t is realized, disordering the sample will 
produce a much more pronounced effect on the in- 
terplane thermal conductivity than that suggested 
by a hybrid gap. 



X. CONCLUSIONS 

In conclusion, we investigated the constraints placed 
on the shape of the superconducting gap in Celrlns by 
the results of Shakeripour et al2i, assuming a one-band 
model with the Fermi surface open along the c-axis of the 
crystal. We find that the temperature evolution of the 
anisotropy between the transport along the c-axis and in 
the plane is incompatible with a simple d-wave gap over 
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FIG. 12: Increase of the residual thermal conductivity with 
impurity scattering rate. 



a quasi-two dimensional Fermi surface, such as often as- 
sumed in the studies of related compound, CeCoIn.5. The 
temperature variation of the anisotropy shows that the 
c-axis conductivity is affected more significantly than its 
intraplanc counterpart by the opening of the supercon- 
ducting gap. This implies that the ability of the unpaired 
quasiparticles to carry heat current along the c-axis is 
impaired relative to their contribution to the in-plane 
transport. This disparity can arise either from the gap 
anisotropy or the Fermi surface properties. We explored 
several models where such modulation has different phys- 
ical origin. 

We find that a weakly c-axis modulated d-wave gap 
does not easily provide satisfactory agreement with ex- 
periment, although may do so with fine tuning of the 
parameters. Of the models with constant interplane hop- 
ping, that with horizontal line nodes most naturally 
account for the anisotropy, in partial agreement with the 
argument of Shakeripour et al. 

Remarkably, we also find that an alternative model, 
where the c-axis hopping depends on the direction of 
the in-plane quasiparticlc momentum, yield results that 
agree with the experimental data even if purely d-wave 
gap, with vertical line nodes, is assumed. Such an 
agreement requires substantial variation of the hopping 
between the nodal and the antinodal directions, in qual- 
itative agreement with the analysis of the Fermi surface 
obtained in the dHvA measurements. 

Given that there is some evidence for a distinct su- 
perconducting dome in Celrlns it is clearly very impor- 
tant to determine the shape of the gap. We therefore 
suggested several possible experiments aimed directly at 
distinguishing the two situations and hope that further 
analysis will help resolve the uncertainty regarding the 
gap symmetry in this material. We also believe that the 
analysis above is generally relevant to a number of uncon- 
ventional superconductors with a quasi-two-dimensional 
parts of the Fermi surface. 
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